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ABSTRACT 

We  report  on  ruby-laser  scattering  measurements  of  Langmuir 
turbulence  excited  in  a high-density  plasma  by  a CO^  laser 
pump.  The  observations  can  be  explained  by  a new  treatment 
of  the  convective  electron-ion  decay  instability  at  equal 
temperatures. 
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Recently,  there  has  been  a great  deal  of  experimental 
1 2 

activity  * devoted  to  plasma  properties  in  the  vicinity  of  the 

critical  surface  of  laser- Irradiated  high  density  plasmas.  Among 

the  diagnostics  that  have  been  employed  are  single-laser  reflection, 

2 

and  harmonic  light  emissions.  In  this  letter  we  report  direct  ob- 
servations of  a decay  parametric  instability^ driven  by  a TEA  CO 
laser  pump  near  the  critical  surface  of  a separately  created  high 
density  plasma.  The  principal  diagnostic  is  Incoherent  scattering 
performed  with  a second  (ruby)  laser.  We  believe  that  incoherent 
saatterlng  can  offer  added  flexibility  and  control  in  comparison 
with  other  diagnostic  methods  for  measuring  Langmuir  turbulence  of 
laser  origin  in  high  density  plasmas. 

In  addition,  a theory  of  the  ruby-scattering  cross  section 

is  set  forth,  based  on  the  linear  convective  electron- ion  decay- 

3 4 

instability  at  equal  temperatures.  * The  theory  appears  to  account 
for  most  of  the  behavior  of  the  cross  section  as  a function  of  CC>2 


intensity. 
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The  deuterium  plasma  is  created  In  a Mather-type  dense  plasma 

focus  device.^  The  measurements  are  made  after  the  current  sheet  has 

run  down  the  co-axlal  electrodes  and  the  current  has  become  mainly  axial 

near  the  anode. ^ In  this  "pinch"  phase,  the  plasma  resembles  a Z-plnch, 

with  density  decreasing  radially  from  a maximum  at  r - 0.  Earlier 
4 

experiments  in  this  device  have  tended  to  show  that  the  plasma  is 
essentially  quiescent  In  the  radial  direction. 

The  laser-plasma  geometry  is  shown  in  Fig.  1.  In  Fig.  1(a)  it 
is  seen  that  the  CO^  and  ruby  beams  are  arranged  co-axially  with  each 
other  (AB  axis),  and  orthogonal  to  the  pinch  axis  (CD).  The  peak  CO^ 
output  of  4.5  J average  delivered  in  40  nsec,  is  focused  by  a 27  cm  Ge  lens 
to  a spot  size  of  750  ym.  The  5J  output  of  the  ruby  beam  is  delivered 
in  15  nsec,  and  focused  with  a 75  cm  lens  to  a spot  size  of  0.25  cm. 

The  collected  scattered  radiation  lies  in  an  annulus  with  scattering 
angles  between  5.5°  and  10.  ^ The  mask  selects  out  only  the  radially- 
propagating  Langmuir  waves  normal  to  the  plane  ABCD.  This  is  also  the 
direction  of  plane  polarization  of  the  CO^  beam. 
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The  scattered  ruby  light  passes  through  an  interference  filter 
centered  at  6500  X with  50  X passband.  If  the  Langmuir  frequency  were 
exactly  equal  to  the  CO^  frequency,  the  plasma  line  would  be  centered 
at  6515  X.  The  50  X passband  is  too  wide  to  distinguish  between  the  plasma  line 
structure  and  an  offset  from  the  CO^  frequency  by  an  ion-acoustic  fre- 
quency (at  most  a few  Angstroms). 

The  photomultiplier  receives  noise  emission  (mainly  bremmst rah lung) 

from  the  plasma,  as  well  as  scattering  from  the  Langmuir  feature.  In 

the  absence  of  the  C0^  beam  the  Langmuir  feature  is  not  observed.  A 

simple  calculation  of  equilibrium  scattering  from  Langmuir  waves,  in 

the  scattering  volume,  V^q,  associated  with  the  50  X passband,  gives  a 

cross  section  a factor  of  two  below  our  experimental  detection  level.  With 

a sufficiently  intense  CO^  beam,  and  for  the  proper  time-phasing  of  the 

lasers,  the  observed  scattering  can  be  up  to  30  times  the  equilibrium 

level  expected  in  V^q.  We  argue  theoretically  that  the  enhanced  Langmuir 

waves  originate  in  a volume,  V^,  four  orders  of  magnitude  smaller  than 

The  CO 2 and  ruby  pulses  are  synchronized  to  peak  at  variou-. 

times  in  the  radial  compression  history  of  the  plasma.  The  Langmuir 

waves  are  only  observed  about  35  nsecs  prior  to  peak  compression 

(-35  nsecs).  At  this  time,  the  maximum  electron  density  is  equal 
19  -3 

to  1.5  x 10  cm  , and  we  assume  the  CO2  radiation  is  reflected.  Radial  and 
axial  density  profiles  are  measured  by  interferometric  techniques,  in 
the  absence  of  the  CO^  beam.  Time-dependent  density  profiles  are 
constructed  over  many  shots  by  Mach-Zehnder  and  holographic  interferometry, 
and  confirmed  over  a single  shot  by  schlieren  photography.  This  tech- 
nique has  been  used  before ^ and  yields  reproduceable  curves  of  n (r,z,t). 
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From  this  we  determine  that  the  reflection  point  of  a normally-incident 

CO^  ray  is  rQ  = 0.04  cm,  and  that  the  electron-density  scale-length 

L = n/|3n/3r|  is  about  0.20  ± 0.05  cm  at  the  time,  -35  nsec,  when  Lang- 
ro 

muir  waves  are  observed.  Geometric  optics  predicts  the  locus  of  CC^ 

ray  reflection  points  indicated  in  Fig.  1(b). 

Figure  2(a)  shows  the  observed  scattering  cross  section  as  a 

function  of  focused  vacuum  CO^  intensity,  at  ”^5  nsec- 

At  times  earlier  than  about  35  nsecs  prior  to  peak  compression, 

the  density  pulse  has  not  yet  converged  to  the  axis,  the  maximum  density 
19  -3 

is  below  10  cm  , and  no  scattering  is  observed,  even  with  maximum  CO^ 

intensity.  During  the  35  nsecs  Just  before  peak  compression,  the  maxi- 

19  -3 

mum  density  is  above  10  cm  , and  L decreases  from  0.20  cm  to  0.025  cm. 

No  enhanced  scattering  is  observed  at  peak  compression. 

At  peak  compression,  the  electron  temperature  is  measured  to  be 

g 

about  2 keV,  by  light  scattering  and  by  spectral  line  emission  . At 

-35  nsecs,  the  electron  temperature  is  calculated  by  assuming  two- 

dimensional  adiabatic  compression.  The  local  electron  temperature  is 

870  - 200  eV.  From  this  and  the  mean  scattering  angle,  0g  = 8°,  we  calculate 

the  scattering  parameter  a = k^/k  - 11. 

9 -2 

The  experimentally  observed  CX^  "threshold"  of  5 x 10  W cm 

is  about  an  order  of  magnitude  above  the  threshold  for  the  homogeneous 

electron-ion  decay  instability  in  the  presence  of  equal  electron  and 
3 4 

ion  temperatures.  ’ We  show  next  that  a linear  convective  saturation 
theory  of  this  instability  is  adequate  to  account  for  most  of  the 
experimental  observations.  The  appropriate  geometry  is  similar  to  that 

Q 

considered  by  Perkins  and  Flick,  but  their  theory  is  only  applicable 


b 


when  >>  T^.  In  our  VTKB  theory  we  employ  a kinetic  model  of  ions  and 

4 

take  into  account  the  important  "beat"  spontaneous  emission  of  Langmuir 

waves.  The  details  of  the  theory,  and  arguments  for  the  absence  of 

10 

other  Langmuir  instabilities  are  given  elsewhere. 

The  integrated  scattering  cross  section  is  given  approximately  by 

a - (r  2/a2)n  V I , V = 2A  , (1) 

e oss  s s 

2 -26  2 

where  r « 5 x 10  cm  . The  solid  angle  for  scattering  is  = 

_3 

6.5  x 10  steradians.  n is  the  mean  local  electron  density.  V is 

o s 

the  scattering  volume,  and  l is  the  axial  length,  equal  to  the  CO^ 

spot  size,  750  pm.  The  cross  sectional  scattering  area  A is  shown 

s 

in  Fig.  1(b).  I( k)  = lim< | 6E(k) | 2>/4ttV0,  is  the  spectral  intensity 

V-K) 

of  longitudinal  field  fluctuations  around  the  plasma  frequency,  normalized 
to  be  one  in  equilibrium. 

For  theoretical  simplicity,  we  straighten  out  the  curved  scat- 
tering volume  in  Fig.  1(b)  into  the  slab  geometry  shown  in  Fig.  2(b). 

The  x-direction  is  the  direction  of  an  assumed  linear  density  gradient 
with  scale  length  L.  The  C02  field  is  in  the  y-direction.  In  the 

Langmuir  ray  paths  in  Fig.  2(b)  k^  is  assumed  to  be  a constant,  set 
by  the  scattering  geometry  to  be  k_/a.  Also,  k - -w  t/2L,  and 

D X O 

x - x ■ -(3/4)v  2t2/L.  u)  is 
o e o 
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the  CO^  frequency,  and  xq  is  the  Langmuir  wave  reflection-point. 

Magnetic  field  effects  can  be  ignored. 

Along  a ray  path,  the  steady-state  kinetic  equation  for 

I “ I (k  (t) , x(t);  k , to  ) is  given  by10 
x y l 


(d/dt)Z  » 2 r / + (// 1)  + 2y  + (w  /w  -<i>  )2y  . (2) 

g c p o L g 


The  first  term  on  the  right-side  is  the  net  growth  rate,  T = v -y  . 

g g c’ 

3 

where  is  the  appropriate  parametric  growth  rate,  and  Yc  is  the 
collisional  damping  rate.  The  second  term  corresponds  to  first  order 
WKB  swelling  of  Langmuir  waves  due  to  the  divergence  of  the  group 
velocity.  The  third  term  is  linear  emission,  and  the  fourth  term  is 
the  dominant  "beat"  spontaneous  emission,  in  which  low  frequency  ion 
Cerenkov  emission  is  upshifted  in  frequency  by  the  CO^-pump,  and  be- 
comes a source  for  Langmuir  waves. 

We  integrate  (2)  from  an  initial  time,  t . , at  which  T =0 

1 g 

(with  the  boundary  condition  /( t^)  = 1),  to  a final  time,  t^: 


Z(tf)  - (tb/tf)[a(mi/me)1/2/1.7]exp[(5/256)Lk  (E^/nG)]  . (3) 


The  exponentiation  factor  depends  on  the  local  mean  (swollen) 

_ 2 - 2 
C0_  field,  E , which  is  related  to  the  vacuum  intensity,  by  E = 

^ o o 

(16ira//3c)IVAC.  Exponentiation  has  been  allowed  to  continue  to  the 

1/2 

reflection  point.  The  factor  a(m./m  ) / 1 . 7 arises  from  beat-emission. 

i e 

The  factor  (t^/t^)  is  due  to  Langmuir  wave-swelling,  which  is 
allowed  to  continue  up  to  the  time  at  which  WKB  breaks  down, 


8 


t a 1 i/u  * “2(k  L)~  /u  /3.  t is  the  effective  time  at  which  swelling 
f WKB  D o b 

10 

begins  in  the  beat-emission  source  term.  Numerical  integrations  give 

9 1 / / 

t.  on  the  ordet  of  t /2»-(Lk  /u>  ) (E  u)  / 32vn&y  ) , and  t,/t,  between 

b i yooo  c bt 

3 and  6. 

The  scattering  area  can  be  shown*®  to  be  A * 266  r 6r,  where 

s s o 

60^  is  the  range  of  accepted  scattering  angles  (5°).  6r  is  determined 

by  the  time  tj^  at  which  J(tj^)  = I(t^)/2,  due  to  the  combined 

effects  of  diminished  swelling  and  exponentiation.  Numerical  results 

show  that,  for  our  parameters,  ~ B'  The  ray  e9uation  yields 

-1  1/3 

for  this  time,  6r  = 6x  = 2k^  (k^L)  > which  is  numerically  on  the 

order  of  5 x 10  4 cm. 

Evaluating  all  the  quantities  in  Eq.  (3)  with  the  experimentally 
determined  scale  length  of  L = 0.20  cm  ± 0.05  cm,  a temperature  of 

18  — 9 

800  eV  ± 200  eV,  and  the  frequency-matched  density  of  9.7  x 10  cm  , 

gives  the  theoretical  curves  shown  in  Fig.  2(a). 

The  nonlinear  saturation  time,^  tM  , due  to  induced-scattering 

N L 

off  ions  in  a homogeneous  plasma,  can  be  estimated  from  the  time- 

12 

asymptotic  analytic  solution  and  compared  with  the  time  available 

along  the  ray  path.  We  find  the  two  times  are  both  equal  to  0.04  nsocs 
10  -2 

when  I =10  W cm  Thus,  it  is  reasonable  that  the  curve  of 

a(IVAc)  should  cease  to  be  exponential  much  beyond  this  intensity,  and 

10  -2 

we  anticipate  that  the  point  at  2 x 10  W cm  requires  a nonlinear 
inhomogeneous  theory.^ 
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The  calculated  threshold  and  its  dependence  on  scale  length  is 

consistent  with  the  experiment,  and  the  inhomogeneous  linear  saturation 

theory  appears  to  give  a satisfactory  fit  to  the  data  at  CC>2  intensities, 
10  -2 

lyAc’  around  10  W cm  and  below. 

One  of  us  (MVG)  would  like  to  thank  D.  F.  DuBois,  D.  Forslund, 
and  D.  Nicholson  for  important  conversations. 
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FIGURE  CAPTIONS 

Fig.  1.  (a)  Scattering  geometry;  (b)  enhanced  scattering  volume. 

Fig.  2.  (a)  The  ruby  scattering  cross  section  as  a function  of 

the  vacuum  focused  CO2  intensity.  Both  experimental 

points  and  theoretical  curves  are  shown.  The  solid 

curve,  C,  is  for  L = 0.20  cm,  and  T * 800  eV.  The 

e 

dashed  curve  C+  (C  ) results  from  an  increase  (decrease) 
of  L by  0.05  cm  or  a decrease  (increase)  of  Tg  by  200  eV. 
(b)  Langmuir  ray  paths. 
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